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Abstract 

In this set of lecture notes we review the mode-coupling theory of the 
glass transition from several perspectives. First, we derive mode-coupling 
equations for the description of density fluctuations from microscopic con- 
siderations with the use the Mori-Zwanzig projection operator technique. 
We also derive schematic mode-coupling equations of a similar form from a 
field-theoretic perspective. We review the successes and failures of mode- 
coupling theory, and discuss recent advances in the applications of the 
theory. 



1 Important Phenomenology for MCT 

Since our objective will be to sketch a derivation of what we will call Mode- 
Coupling Theory (MCT), we will focus our attention on one observable in 
particular, namely density fluctuations. For this we will first define some 
of the concepts needed to do so. 

We want to calculate a specific time correlation function. In general, 
such a function is expressed as follows, 

c{t) = {Amm- (1) 

It is an ensemble average of the evolution of the fluctuations of a variable 
in time, at equilibrium. As seen in Fig. A{t) fluctuates around its 
average value in equilibrium, while C{t) measures the correlation of A at 
one time with the value of A at another time. 

The density or particles in a liquid can be one example of A{t), 



,(r,t) = ^5(r-r,W), (2) 



* Current address: Department of Chemistry, Columbia University 3000 Broadway, New 
York, NY 10027 reichmanOchem. columbia.edu 
tpcharbon@fas.harvard.edu 



1 




L 



Figure 1: Left: time evolution of the instantaneous fluctuations of the quantity 
A. The product of fluctuations separated by time t' , averaged over all ^o's gives 
the correlation function C{t'), at equilibrium. Right: fluctuations of the density 
on lengthscale L ^ 27r/fc; if k is small, the area of density fluctuations is large. 

which we can Fourier transform, 

Pk(t) = /rfre* '"<5(r-r.(t)), 

i 

= ^ (3) 

i 

In this case the correlation function will be labelled _F(k, t), which is can 
be expressed as follows. 

Fik,t) = ^ (p-k(O)pkW) = ^ E (e-'-'fo'e*-^'*)} . (4) 

ij 

Note that we need to have X^i ki = (i.e. — k + k = 0!) to conserve 
momentum, otherwise the correlation function is equal to zero. 

The variables labeled by k measure density fluctuations in reciprocal 
("fc = |k|") space, which can be thought as the inverse length. When k is 
small we are looking at long lengthscales, as we can see in Fig. Q When 
it is large, we are probing very short scales. 

The function F{k, t) is essentially what scattering experiments mea- 
sure. At t = 0, 

F{k,t = 0) = 1 {p-k(0)pk(0)) = S{k), (5) 

where S{k) is called the static structure factor of the liquid. Why that 
name? Consider the radial distribution function of a liquid g{r). The 
function g{r) is proportional to the probability that a particle is a distance 
r away from a particle at the origin. In a dense liquid g{r) shows the 
structure of the solvation shells as depicted in Fig. |5| Also, it can be 
shown that P H E| 

S(fc) = l + p /dre-*-<?(r), (6) 
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Figure 2: Left: radial distribution function g{r) for a simple liquid of size a. 
Right: the corresponding structure factor S{k). A sample structure is also de- 
picted where the solvation shells are indicated by the dotted lines. The exclusion 
radius can be seen in the absence of amplitude of g{r) for r <C cr. 

where p = N/V is the density of the system, and thus S{k) is also indi- 
cating something about the hquid structure 0121 01 • An example for a 
simple liquid is depicted in Fig. |21 

But how do we expect F{k, t) to behave? For high temperatures - 
above the melting point - F{k, t) will decay like a single exponential func- 
tion in time for for k > 2TT/sigma as plotted in Fig. (31 For supercooled 
liquids, the situation is different and a characteristic decay pattern can 
also be seen in Fig. El We observe a multi-step relaxation. 

1. At short times decay is coming from free and collisional events that 
involve local particle motion. Consistent with a short-time expan- 
sion, F{k, t) ~ S{k) - A{k)t^ + ... in this regime 0(51121 • This will 
be true at any temperature. We will not be concerned much with 
this part of the decay. 

2. Intermediate times encompass a period during which particles ap- 
pear trapped in cages formed by other particles. This regime is the 
/9-relaxation regime. The decay to the plateau (Ila) may be fitted 
as / + At~'^ and the decay from the plateau (lib) as / — Bt'^ . Also, 
the exponents have a scaling consistent with the relationship 



r(l -g)^ _ r(l + b)^ 
r(l-2a) " r(l + 26)' 



(7) 



3. At long times, in the a-relaxation regime, the decay may be 
to a stretched exponential law 4 



fitted 



F{k,t) ~ e-^^y . 



(8) 



with < /? < 1. Do not be confused with the notation. It is the /3 
power that appears in the a-relaxation regime! In general {3 and r 
will be k and temperature dependent. 
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Figure 3: Left: F{k,t) exhibiting exponential e~*/^ decay for a normal liq- 
uid. Right: supercooled liquids do not have such a simple decay. The various 
temporal regimes are described in the text. Notice the logarithmic scale. 



In a later section, we will return to this kind of phenomenology. For 
now, we just make a few superficial remarks about things that will be 
covered in more depth by others in these lectures. 

The constant r that appears in the stretched-exponential decay law 
is strongly temperature dependent. All transports coefficients -D (diffu- 
sion), rj (viscosity), etc. - are strongly temperature dependent as well. 
Over a rather wide range of temperatures, a fit to this temperature de- 
pendence may be 0] 

^^e'^. (9) 

Clearly, as To is approached, relaxation times become so large that the 
system cannot stay in equilibrium. Other fitting forms, some that do not 
imply a divergence at finite temperatures, may be used to fit the data as 
well. 

Some systems, hard-spheres for example, are not characterized by tem- 
perature, but by density of packing fraction = ^na^p, where a is the 
particle radius. For such systems, one may fit with |1| 

rj , (10) 

or with other forms. 



2 The Mode- Coupling Theory of density 
fluctuations 

Our strategy will be to derive an exact equation of motion for F{k, t) 
and then to make approximations that allow us to solve them [5)1^1. The 
approximations are uncontrolled, and we will judge them by their success 
or failure. 
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2.1 Memory functions 

Consider some classical function of phase space variables A{t), where the 
time dependence originates from that of the positions and of the mo- 
menta Pi for a Ai'-particle system. We know from Hamilton's equations 
that 

^ = {A(t),H}^ii:A{t), (11) 
where {, } is a classical Poisson bracket, which can be expressed as follows, 
r . oi V- fdA dB dA dB\ 

Also, for liquids of interest, 7i is a classical Hamiltonian with pairwise 
interactions ^(r) between the particles, 

^ - E + ^ E ^(r^.)- (13) 
We can thus identify the following. 

It would be possible to integrate the differential equation to find A{t) = 
e^''*A(0), but this is not useful by itself. 

We also need to define a scalar product of the variables as 

iA,B) = {A*B). (15) 

Now, consider an operator called a projection operatorm^] |21 130 

V = {A,...){A,A)-'A. (16) 

If j4 is a vector, {A,A)~^ is thus the inverse of a matrix. Note also that 
P^A — VA = A. In geometrical terms, the projection operator finds the 
component of some variable B along the chosen direction A, as depicted 
in Fig. H 

This is useful as we can extract from an arbitrary B how much "char- 
acter" of A it has. In particular, the operator A may be a slowly varying 
(quasi-hydrodynamic) variable. Consider the density as defined in Eq. 

i 

and then 

Pk(i) = ik . ^ PlWe*-r.(t) ^ ,k . hit) = i|k|j^ (t), (17) 

i 

where jk(^) is the longitudinal current. If k is small (large lengthscales) , 
then pu.{t) is approximately small. This is what is meant by slow. In the 



5 




Figure 4: A two-dimensional version of the projection operator VB. The quan- 
tity B is projected unto the space A, which extracts the A component of B 
(indicated by a thick dashed hne). 



limit A; = 0, then pk^o = 0, and the density is strictly conserved. As Fig.0 
indicates, if k is small, the area of density fluctuations if large, i.e. the 
rate at which the number of particles fluctuates is small. 

We now want to find the exact equation of motion for a correlation 
function 

^ = e'^''[r + ii~r)]iCA, 

= in-A{t) + e'^\l~P)i£A, (18) 

where 



ig,= {A,iCA) ■ (A,A)-\ (f9) 



Now, writing 



e — e 



0{t) can be obtained by differentiating both sides of the equation 



iC ( e'^^^Oit) + e 



iLe 




i[{i-r)c] e*(i- 


-V)Ct 

1 


i(\-V)Ct\^ 




j[(f-P)£] e*(i- 


-V)Ct 

1 




= e'^'e)(f) + j[(f-P)£]e 


i(l-V)Ct 












0{t) 


rt 

■ 1 7 ijCr^-rt i(l — 'P)CT 

— I are VCe ^ ' 
Jo 
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(21) 



where the last equality follows from the fact 0(0) =0. As a result, we 
may write 

e''^\{\-V)LA = /"dre'^('-"'jp£/(r) + /(f), (22) 

where /(t) is called the fluctuating force. 

/(f) = e'<'-^)^'j(f-P)£A (23) 
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What does this mean? The fluctuating force is obtained by taking the 
time derivative of A, using the complimentary projection operator (1 — V) 
to remove the "A" character - perhaps the slow character - and is then 
propagated in the orthogonal - fast - space. To put it another way, if V 
removes the slow character from a variable, then the fluctuating force is 
the remaining fast force. We will come back to this later. 

It can be shown that [A, f{t)) — 0, by noticing that the deflnition of 
f{t) contains the (1 — V) factor. This means that f{t) is orthogonal to A, 
in accord with the discussion above. Noting that 

iiA, £f{t)) = i{£A, fit)) = - P)£A, fit)) = -(/(O), fit)), (24) 

the first term on the RHS of Eq. 1221 allows the equation of motion to be 
rewritten as 



dAit) 
dt 



iQ ■ Ait) ■ 



drKir)-Ait-r) + fit), 



where we define the memory function Kit) 

Kit)^if,fit))-iA,A)~\ 



(25) 



(26) 



This is a fundamental and exact equation for the time dependence of Ait). 
Defining the correlation matrix 



Cit) ^ (A* io) Ait)) = iA, Ait)) 
and using the equality (A, fit)) — 0, we get 
dCit) 



dt 



= in • Cit) 



drKir) ■ Cit ■ 



(27) 



(28) 



as an exact equation for the matrix of correlation functions C(t) that we 
will want to compute. The problem with computing Cit) is embodied in 
the difficulty of determining Kit). 

Now, we want to focus this general framework on density fluctuations 
with respect to the bulk density p, which will allow us to get an expression 
for the intermediate scattering function. Consider |3 



A = 



(29) 



where 



■ L 



i 



(27r)V5(q), 



(30) 



and therefore 



Cit) = 



(A'Ait)} 

(Sp-qSpqit)) 

U-^Sp^it)) 



{Sp-^Ai^)) 



(31) 



7 



For the purpose of this demonstration, we will concentrate on the clement 
in the lower left corner of the matrix, which is in this case ^ dt^'*^ • -^t 
t=0, the matrix reduces to 



00) = 



NS{q) 

Q NkBT 



(32) 



Also, 



in = {A* A) ■ {A*A)-\ 





;NqkBT_ 

m 



.■ qksT 

'mS(q) 



c 







JVS(q) 









NkBT 



(33) 



To obtain the last two equations, we used integration by parts, the prop- 
erty that the correlation of an observable and its derivative is always zero, 
and the statistical thermodynamics result that follows: 



(jVp,) = -^^((q.p,)e-^<'-*(q.p,)e^''-^), 

.^{jnv,) =1——. (34) 



iq 
m 



The random force /(O) is expressed as 
/(O) = {1-V)A, 





■ qkBT_ 
' mS{q) 



dt 





' " 







We will now look at the equation of motion term by term. First, 
Note that the lower left corner term equals -5^ rffi?'*^ • Second, 



iQ ■ C{t) = 





i qksT 

mS(q) 



(<5p-q5pq(t)) (6p-^j^it)} 



(35) 



(36) 



(37) 



Notice the lower left corner term is —^^^^F{q,t) this time. Lastly, the 
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memory matrix is 



Kit) = 



[0 





{R-qRq{t)) 



NS{q) 









(38) 



Concentrating on the lower left corner, using the equation of motion from 
Eq.|2Hl we find PH|H1|S1 



<fF{q,t) q^ksT 



+ 



NknT 



j\r{R.Mr))^/{q,t~T)^0. (39) 



This equation is exact, but impossible to solve. To make approxima- 
tions, we will look at {_R_q7?q(f)) using some intuition. Recall that 



dt 



. qksT 
mb(q) 



(40) 



and that 



dt 



d 
dt 



1 v-^ /- rfP 



dt 



' + ;^E(^-pO'^'"^'- (41) 



Also, note that in this last equation is a force and therefore 
^ ~ ^E ^''^(l'"' ~ '■^D ^ ^ ik(?!)k5pk(5p-k. 



(42) 



where we made the momentum-space transformation 4>\^ — j dre^^"^ 4>{r). 
We then see that hidden in the fluctuating force is a pair of densities. This 
illustrates an important point: at first, we may have suspected that the 
fluctuating force is a fast variable and that it decays on a short timescale 
because we removed the slow modes 5pk from it - but we see that it 
contains, at leading order, a slow character (at least if 5pw is slow) from 
the product of slow modes 5pk(5/9-k! Overall the time derivative of the 
current has the symmetry of 5pk<5pq-k, where the q factor comes from 
e"' ''' that multiplies the force in Eg.ETI 

We will now approximate (i?_qe'®''"*i?q). As a convention, note that 
Q = l-V. 

-p2e''^*p2, where we define the new projection op- 



1. Replace e 
erator 

V2 



iQCt 



E ^ki,k2(^k3,k4 • • •)(^ki,k,^k3,k4)"' (43) 



ki,k2,k3,k4 
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and where j4ki,k2 ~ 5p]f^^5p\i^. The V2 operator simply projects 
7?q onto its dominant slow product mode. We neglect the Q in the 
exponent, simply because it is hard to compute anything keeping it 
there. However, to 0{q^), it may be shown that this neglect is not 
consequential. 

2. Factorize four-point density terms into products of two-point ones. 
Using this algorithm, we get 

P2i?q= Vq(ki,k2)<5pki<5pk2, (44) 

ki,k2 

where 

V"q(ki,k2) = ^ ((Jpki'Jpka^q) ■ ('^Pki'Jpka'Jpka^Pki)"^ (45) 
k3,k4 

The denominator has a product of four density variables that will be 
factorized into products of two structure factors. The numerator has 
terms like: 

5p__k(5pk-q^^ = -((5/9_k<5pk-qiq) " ((5p-k(5/9k-^q jq ) , (46) 

where the result was obtained by integration by parts, and 
Let's calculate one of the terms in Eg. 1461 

- (5p_k5pk-qJ^) = z/^(k.pOe-*--^e'(''-'')-'^(q.pOe'<^-' 



ij 

= i(k.q):^iVS(|k-q|), (48) 

where we used the result of Eg. 1341 to complete the calculation. The other 
term similarly gives 

- {5p-k<5pk-qJ^> = i (q ■ (q - k)) t^NSik). (49) 

The term in Eg. l47l is hard to compute directly, but within the convolution 
approximation, it can be reduced as follows 01^1: 

(c5p-k5pk-q5pq) ~ NSik)S{q)S{\k - q|). (50) 

After treating all static density fluctuations within the Gaussian (and con- 
volution) approximations, we find that the vertex 'Ki(ki,k2) can be ex- 
pressed as a function of only two wavevectors. As a consequence of transla- 
tional invariance, we would be left only with terms involving the difference 
of wavevectors k = ki — k2, which allows us to write V!q(ki, k2) = Vk,q-k- 
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Also note that the summation is now only over k. Combining all terms 
gives 

_ ikBT \ (q-k) q-(q-k) \ 

= ^5^{(q-k)c(fc)+q.(q-k)c(|k-q|)}, (51) 

where we have rewritten the result using the direct correlation function 
c(fc) = i ^1 — • So, piecing this together, 

((i?qP2)*(P2i?q(t))) ^ |^k,q-kl4',q-k'| (<5p„ k' ^Pk' -q ^Pk (f ) <5pq- k (t ) ) , 

k.k' 

^ El^k.q-kVk'.q-k'|^'nfc,*)mq-k|,t)(5k,k' +4'-q,k), 

k,k' 

= ^'amP' E |v-q-k,k|V(fc,t)F(|q-kU), (52) 

k 

where we used Wick's factorization and where 

l/q_k,k = {(q • k)c(fc) + q ■ (q - k)c(|q - k|)} . (53) 

We need to convert the discrete sum to the continuous integral, ~» 
J dk, and multiply by the prefactor as obtained in Eq. 1391 to 

get the final MCT equation. 



= 


dt^ mS(q) Jo 


,dF{q,T) 
dr 


(54) 


with 








K{q,i) = 


[ dk|K,_k,k|V(fc,t)7^(|k-q|,t) 




(55) 



2.2 Some Properties of the Solution(s) of the MCT 
Equation 

2.2.1 Schematic MCT 

Via the approximation discussed by Bengtzelius et al. we can reduce our 
MCT equation to a schematic form, 

^ + mm + X /J dr<,\t -rf-^^O, (56) 

where ^{t) ^ F{k,t), before we neglect the coupling wavevectors. The 
solutions of Ea. l56l have been discussed by Leutheusser \7\ and Bengtzelius 
et al. (2- The most striking feature of this equation, and of the full MCT 
equation from which we "derived" it, is that there is a transition to a 
completely non-ergodic phase for particular SIq ^md A (or T and p for the 
real MCT equation). The two cases are depicted in Fig. El 

The transition from ergodic to non-ergodic at a sharp, well-defined set 
of parameters may be interpreted as the transition from a liquid to a solid. 
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Figure 5: Left: decay in the ergodic (supercooled) case. The correlation 

vanishes on a finite timescale. Right: in the non-ergodic (glassy) case, that same 
function remains finite even for infinite times. 



The fact that correlations do not decay as t ^ oo is indicative of this. 
However, no information of an ordered state was used or imposed. Thus, 
the solid could only be a disordered one, i.e. a glass. 

2.2.2 Solutions of Full MCT Equations 

If we denote by Tc the temperature where MCT predicts a glass transition, 
the relaxation times r scales as |S||^|7||3 

T{q,T)^A,{T-nr'', (57) 

which means that when T Tc, t diverges as a power law with a universal 
exponent 7. This form may fit data, but only over a limited temperature 
range. 

Also, the decay in the /3-relaxation regime is indeed given by [Sl|ni|71|3| 

early 13: f + At-"", 

late /3: / - Bt\ 

with 

r(l~2a) r(l + 26)' ^ ' 

at least for T very close to Tc, i.e. when ^ 1. This is a great 

triumph of the MCT equations and is fully consistent with simulations and 
experiments Furthermore, the power 7 is related to the exponents 

a and b as 

For the a-relaxation regime, an approximate solution of the full MCT 
equations is indeed approximately given by 

F{k,t) ~ e V-J , (60) 
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and this stretched exponential function describes well experiments and 
simulations. The schematic equation, however, only exhibits exponential 
decay. 

More generally, MCT predicts that for a correlator at temperature T 
a time-temperature superposition holds 

C{t,T)=Cit/T{T)), (61) 

where C{t, T) is a correlation function, C is some master function and t{T) 
is the a-relaxation time. This is also generally consistent with experiments 
and simulations. 

2.2.3 Redux: An Assessment of the Successes and Failures 
of MCT 

Even a scientist who is opposed to the spirit and approximations that go 
into MCT ought to be impressed by its success, where it succeeds. Fur- 
thermore, it is essentially the only first-principle theory of glassy liquids. 
Namely, from the structure of the liquid alone {S{k), the structure factor) 
a detailed set of dynamical predictions emerge. We will now spell out 
where MCT works and where (we think) it does not. 

Successes 

1. MCT makes some remarkable predictions that are correct. For ex- 
ample, the remarkable scaling properties in the /3-relaxation regime 
that are predicted are essentially correct and so is the time-temperature 
superposition in the a-relaxation regime. This is similarly accurate 
for other predictions that we will not discuss here 5 6 . 

2. MCT has predicted novel relaxation patterns correctly. One recent 
striking example is the behavior of colloidal suspensions with induced 
short-ranged attractions|H|. Here, MCT has predicted that adding 
attractions may melt the glass from hard-spheres, and that for cer- 
tain parameters, logarithmic relaxation may be observed. Both pre- 
dictions have been confirmed again by computer simulations and 
experiments. 

3. Although we will not discuss it here, there exist models with quenched 
disorder (spin-glass models) for which the schematic MCT is exact 
|1UI lll| . These models make connections between MCT and energy 
landscape theories possible, as well as extensions of the MCT ap- 
proach to situations that are out-of-equilibrium (aging). This has 
been very fruitful and has led to new insights into glassy systems. 

Failures 

1. The best-known failure of MCT is that it predicts a sharp glass 
transition at a temperature Tc, but Tc > Tg. This means that MCT 
predicts kinetic arrest to a non-ergodic phase at temperatures where 
the system is still ergodic and liquid. 
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Figure 6: The non-Gaussian parameter as depicted here is one of the features 
that MCT cannot reproduce accurately. The curves increase in magnitude and 
spread with decreasing temperature (from left to right). 



2. MCT also predicts power-law divergence of transport coefficients 
and relaxation times as in Eq. 1571 but this is less accurate over 
a wide range of temperatures than the temperature dependence of 
transport coefficients given in Eq. |5| It is a reasonable fitting 
form over several decades of relaxation time in mildly supercooled 
liquids. In addition, the parameters (3, a, b, . . . are predicted to be 
constant in MCT, at least for temperatures for <C 1. But in 
actuality, they are mildly temperature dependent. One should be 
careful, however, not to take asymptotic predictions of MCT and 
apply them to cases where ig not small 5 . 

3. Another failure of MCT is in the prediction of certain indicators of 
collective relaxation. In general, timescales and lengthscales of such 
heterogeneous motion can be probed by multi-point correlations. A 
simple, non-multi-point function that seems to correlate crudely with 
the timescale of such motion is the non-Gaussian parameter 

where 

{r\t))^{[r{t)~T{Q)f) (63) 

for a tagged particle and similarly for the other term. Usually, the 
behavior of {a^{t)) is similar to what is depicted in Fig.|^ However, 
MCT predictions of a2 {t) are quite inaccurate. In general, MCT fails 
to accurately depict properties such as the non-Gaussian parameter 
and the breakdown of the Stokes-Einstein relation. 

One may take from this last result that MCT is not capable of saying 
anything about dynamically heterogeneous motion in supercooled liquids 
|16l . but perhaps this statement is too strong. We will explore this further 
in the next section. 
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2.3 Field-Theoretic Description 

Before finishing this section, we provide a sketch of the field-theoretic 
approach to schematic MCT |1UII11| . In some sense, the memory function 
approach can be thought of as arising from coupled Langevin equations 
for the modes iSpq and jq . For example, 

= ^gjqW, (64) 



from the fluctuating force 



- ^JqW + ^ ■ (65) 

^ ^ ^ noise term 

viscosity term 

As a toy model for this (forgetting vector labels and wavevectors) , we get 

0(t) = -Ki)0(i)"f'^(i)'+^(i) (66) 
with {r]{t)ri{0)) = 2T5it). We also define 

Goit,t') = e-^'"''''^''> (67) 

and 

(G(g)/)(t)= r dt'Go{t,t')f{t'), (68) 
Jo 

where Go is the bare response function, or propagator. It must be zero if 
t < t' . The solution for 0(t), with <^(0) = is, in graphical terms, 

m = (69) 

where an arrow represents Go, an x represents the noise, and a factor of 
g/2 is associated with each branching point. These terms are simply ob- 
tained as a solution from integrating Eq. 1661 In a more compact notation, 

Ht) = Go (g) r, - |Go {Go r; • Go »?} + . . . , (70) 

where "•" is a simple product. 

We define two kinds of functions, 

c{t,t') ^ imHt')) (71) 

and 

where the last equality is true for Gaussian noise. Again, C{t,t') can 

be defined on the entire {t,t') plane but G{t,t') only for t > t' . Let's 
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construct a series for the two functions. The zeroth-order contribution to 
C{t, t') is given as follows: 

= (^J^ dtiGoit,ti)f{ti) dt2Go{t',t2)f(t2) 

= [ [ dtidt2GQ{t,ti)Goit',t2)2TS{ti-t2). (73) 
Jo Jo 

The bracket average implies connecting the x -vertices to form diagrams 
with none of these left. Diagrams that do not pair up all noise vertices 
{i.e. those with an odd number of x-vertices) average to zero. 
Going beyond zeroth-order, we get the following: 

C{t,t') = +() + ... 
= +..., 

= dti dt2Go{t,ti)D{ti,t2)Go{t' ,t2), (74) 
Jo Jo 

where to obtain the last line we defined 

2 

D{ti,t2) = 2TS{ti-t2) + ^Co{U,t2)Co{ti,t2) + . . . (75) 
What about G{t,t')? 

G(t,t')= + (^)+---, (76) 

where the lone x has to be attached to a x on the tree and then the 
remaining diagram is closed. Thus, 

Git,t') = +..., 

= Go{t,t') + J^ dti J^' dt2Go{t,ti)E{ti,t2)Goit2,t'),{77) 

where we similarly defined 

E(ti,t2) = g^Goiti,t2)Go{ti,t2) + ... (78) 

This is an exact, formal representation of the perturbations series. In fact, 
in some sense it is simply a definition of the kernels D and E. We can 
appeal to the structure of the perturbation series to justify this. 

Also, the lower limit of the second integration insures that t2 > t' . In 
fact, let's take a closer at the limits of integration for a sample diagram, 
the second term in Eq. 1771 which is reproduced on the left-hand side of 
Fig. |7| with additional labels. The incoming branch imposes t > t\, the 
central loop t\ > t2, and the outgoing branch t2 > t' , for an overall 
t > tl > t2 > t' . The resulting integration limits are thus J^, dt\ J^*/ dt2. 

If you have followed so far, you might be asking yourself what happened 
to the diagrams like the one appearing on the right in Fig. |7| Such closed 
loops are called tadpoles. They do not contribute to the time dependence 
of C or G, and we assume that their contribution is absorbed into /i(t). 
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Figure 7: Left: As an exercise, work out the integration limits for this sample 
diagram. The time pairings indicate the beginning and end times of a given 
segment. There are in this diagram two internal vertices, ti and t2, and two 
external ones, t and t' . Right: a simple tadpole diagram. 

The next lowest-order terms are then the diagrammatic forms we last 
drew in Eg. ITU and 1771 

To make a self-consistent approximation we replace Go and Co in our 
second-order approximations for D and E by G and G. We will call this 
the Mode-Coupling Approximation, a name that will be clear in meaning 
at the end (TnHni- This is Eg. 171 and 1771 with 

T.{ti,t2) = g^C{ti,t2)C(ti,t2), (79) 

2 

D(ti,t2) = 2T5(ti,t2) + ^C{ti,t2f. (80) 
These two equations can be further manipulated by noting that 

Go = (^Kt) + ^) , (81) 

and so we can multiply both sides of Eg . 1741 and 1771 by Cq^ . 

Go'0G = / + S0G, (82) 
Go'(8)G = Go'0{Go + Go0S0G}(g)i30G, 
= D(g)G + E(g){G(g)D(g)G}, 

c 

= Z)0G + E0G, (83) 



where / is the identity operator. In other notation. 



1- 

\dt 


+M(t) 


•G{t,t') 


= 5{t~t') + j dt"T.{t,t")G{t",t'), 


(84) 


I- 

\dt 


+Kt) 


'C{t,t') 


= r dt"D{t,t")c{t',t")+ r dt"E{t,t")c(t",t') 

Jo Jo 


(85) 



This still does not look like the MCT equations derived before. We 
will manipulate the RHS of the second equation to achieve this. Taking 
the first term, we substitute 



l^dt"Dit,t")G{t',t")^ 1^ D{t,t")^^C{t',t"). (86) 

We can replace t ^ t' because of the restriction on time arguments in 
C{t' , t") and we can substitute G for G as written by using the fluctuation- 
dissipation theorem (FDT) and assuming that the system is at equilib- 
rium. The result can now be integrated by parts, 

1 [0'(M")G(i",t')]|o - dt"^D'{tX)C{t",t'). (87) 
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(90) 



Note that we can neglect the 5-function part of D here, since it provides 
no contribution. So, we denote the regular part of D as D' . 
Again, using the FDT, we can substitute 

A^D'{t,t")^mt,t"), (88) 

yielding an integral which can be combined with the other term of Eq. 1851 
to obtain 

- r dt"-E{t,t")C{t",t')+ f dt"T.{t,t")C{t" ,t') ^ r dt"E{t,t")C{t",t' 
Jo Jo Jt' 

Now, using the FDT in the reverse direction and integrating by parts, we 
get 

i [D'{tX)C{t'\t')]\l, ^lj\t"D'{t,t")^C{t",t'). 

Combining all terms on the right-hand side of Eq. 1851 we get 

i [D'{t,t)C{t,t')^D'it,0)CiO,t')]~^ £dt"D'{t,t")^C{t",t'). (91) 

Using the fact that, in equilibrium, these functions are time-translation 
invariant, 

D'{t,t") = D'{t-t"), 

C{t,t') = C{t-t'), 

C{t",t') = Cit"~t'), 

C(0,t') = Cit'), 

D'{t,0) = D'{t), 

D'{t,t) = D'iO), (92) 

and making the transformation, 

t — t = T , 

t-t' = r, 

t' -> oo, (93) 

we find 

+ m) C{t) + drD'it - r)^^ = 0, (94) 

where 

m = Kt)~^D'{0), (95) 

D'it) = ^Cit^rf. (96) 
This is just the schematic model with 

nl ^ m, (97) 

A^-^. (98) 
2T ^ ' 
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The only difference is that ^^p- appears instead of This actu- 

ally makes no difference as far as the glassy properties are concerned. 
In fact, for models of overdamped systems, such as Brownian colloidal 
spheres, ^^7^ is what appears naturally in the reduction of the full MCT 
equations of the schematic model. This completes the relationship be- 
tween the memory function/projection operator MCT derivation and a 
field-theoretic approach. 



3 Looking Ahead: Beyond "Simple" MCT 

In this section we will outline some thoughts on attempts to do better 
than the MCT derived thus far. This is not an exhaustive discussion, but 
is meant to give an idea about what can be done. 



3.1 Coupling to Currents 

Gotze and Sjogren |12| as well as Das and Mazenko |13| have developed 
theories that remove the sharp transition at T = Tq. In both cases, it is 
the coupling to certain current modes that are ignored in the expressions 
we have just derived that restore ergodicity below Tc- 

In the theory of Gotze and Sjogren, the Laplace transform of our exact 
equation of motion for F{k, t) is given as 

nk,z)^—4^, (99) 

^ z + M(k,z) 

where Q,{k) = ^^^^"^ and M(fc, z) is the Laplace transform of the memory 
function 

M{k,z)= dte-''^M{k,t), (100) 
Jo 

and similarly, 

F{k,z)= / dte'''^F(k,t). (101) 



Essentially, within the "extended" MCT of Gotze and Sjogren, 

^(^'-)- i-MMmM) ' ^'''^ 

where K{k, z) is the ordinary MCT memory function given in Eg. 1391 The 
expression for Sh{k,z) is complicated, but to see what it does, note that: 

1. if 5h = 0, we recover exactly the ordinary MCT that we have derived 
before. This can be checked by applying a Laplace transform the old 
expressions. 

2. if 5h{k, z) has no singularities as z ^ Q, the the strict transition at Tc 
in the MCT we have previously derived is removed since relajcation 
is governed by M ~ '^/\Sh\ at long times and not K{k,z), which 
yields a pole singularity in z-space. 
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In the theory of Das and Mazenko, a hydrodynamic approach is used. 
The kinetic energy of the free energy functional in terms of current j and 
density p has the form 

K.E.[j,p]^^-. (103) 
P 

This is hke the usual /2m kinetic energy. However, the 1/p part coupled 
to the current rounds off the strict singularity at Tc, like in the Gotze and 
Sjogren theory. 

Note that in both theories, we need currents to restore ergodicity. For 
some systems, hke (simulated) colloidal hard spheres undergoing Brow- 
nian motions, these currents do not exist! Thus, the Gotze and Sjogren 
and the Das and Mazenko theories cannot be used to improve ordinary 
MCT there. 



3.2 New Closures 

As mentioned above, the extended MCT of Gotze and Sjogren and of 
Das and Mazenko cannot tell us anything (beyond ordinary MCT) for 
Brownian hard-sphere systems. An interesting proposal was recently put 
forward by Szamel |14| . The main idea is not to factorize the memory 
function expression leading to the ordinary MCT given in Eq. 1391 but to 
write an exact equation of motion for it, and then factorize the memory 
function for the new equation. Here is a sketch of the idea. 
Recall our old approach to MCT, 

d^F{k,t) k^ksT dF(k,z 



dt^ mS{k) 



F{k,t)+ [ dTK{k,t~T)^^^^^^ =Q, (104) 
Jo dr 

where essentially K{k, t) ~ (SpSpSpSp) is a four-point function of density 
variables. In the old approach, the closure involved 

K(k,t) ~ {5pSp){SpSp) = ^F(g,t)i^(lk-ql,t) (105) 

q 

and this allowed to solve for F{k,t). 

Instead of factorizing the four-point memory kernel, let's write an exact 
equation of motion for it, following the same lines of reasoning as before, 

d'^K(t) [\ nf .9K{t) 



+rK{t)+ dTR{t^r)—^ =0, (106) 

The wavevector indices are suppressed to simplify the notation in order to 
clarify idea behind the manipulations. This has the same form as before, 
but with new frequencies F and a new memory function R{t — r). 

Schematically, R ~ (SpSpSpSpSpSp) is a six-point function! We can 
close the equation for K and F, by approximating R^i K ■ F the product 
of a four-point and a two-point function. This yields two coupled sets of 
integro-differential equations that may be solved self-consistently yielding 
a converged F{k,t). 

This approach has not been considered for the full dynamics of F{k, t), 
but yields a better estimate for Tc, (i.e. the Tc that is extracted is closer 
to the measured glass transition). 
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Figure 8: The non-Gaussian parameter 0:2 (i) (left) and the four-point correlation 
function X4(t) (right) peak at times t* where t* > i*^. 



3.3 Four-point correlations and dynamical het- 
erogeneities 

It was mentioned in Section |21 that MCT does not describe well 

and that 02 (i) seems to correlate well with the timescale of maximal 
dynamical heterogeneity |16l I17| . It turns out that this timescale is in 
the late /3-regime. This highlights the fact 02 (t) yields information on 
transiently mobile particles that jump due to the destruction of cages. It 
should be noted that there is no lengthscale dependence in a2{t). 

To gain some information about a growing (dynamical) lengthscale, 
a multipoint dynamical generalization of the static structure factor may 
be studied T8 19 20 The — > limit of this structure factor, as 

Sharon Glotzer discusses in her lectures, is the susceptibility 



X4(i) ~ y dri . . . dr46la (|ri - r2|) (|r3 - r4|) X 

(p(ri, 0)p(r2, t)p(r3, 0)p(r4, t)} , (108) 

where the function 6a (|ri — r2!) equals one when |ri — r2| < a, and zero 
otherwise 1181 . The timescale at which X4,{t) peaks is generally in the a- 
regime. The growing lengthscale associated with dynamic heterogeneity 
is associated with slow moving, transiently caged particles. 

Given the superficial similarity with 02 (f), as shown in figur^ one 
might conclude that MCT cannot compute objects like X4, but this is not 
the case. Recent work by Biroli and Bouchaud 17, motivated by the ear- 
lier insight of Franz and Parisi ^20^ and Kirkpatrick and Thirumalai |21l . 
shows that MCT may make quantitative statements about the scales of 
length and time associated with dynamical heterogeneity. So far, absolute 
lengthscales have not been computed, but dynamical exponents z relating 
timescales r and lengthscales ^ have, 

r~r, (109) 
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where z = 27 and 7 is given in Eq. 1591 
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